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■ Abstract 

PL 

ly-j ■ We propose a new graph metric and study its properties. In contrast to the standard 

CN . distance in connected graphs [5], it takes into account all paths between vertices. 

Formally, it is defined as d(i,j) = qu+qjj — qij—qji [11], where qij is the (i,j)-entry of 
the relative forest accessibility matrix Q(e) = (7 + eL) -1 , L is the Laplacian matrix 
of the (weighted) (multi)graph, and e is a positive parameter. By the matrix-forest 
(-h \ theorem, the (i , j)-entry of the relative forest accessibility matrix of a graph provides 

the specific number of spanning rooted forests such that % and j belong to the same 
tree rooted at i. Extremely simple formulas express the modification of the proposed 
distance under the basic graph transformations. We give a topological interpretation 
of d(i,j) in terms of the probability of unsuccessful linking i and j in a model of 

■ random links. The properties of this metric are compared with those of some other 
C*~) . graph metrics [13,1]. An application of this metric is related to clustering procedures 

such as centered partition [2]. In another procedure, the relative forest accessibility 
£\ and the corresponding distance serve to choose the centers of the clusters and to 

assign a cluster to each non-central vertex. Some related geometric representations 
are discussed in [3]. The notion of cumulative weight of connections between two 
vertices is proposed. The reasoning involves a reciprocity principle for weighted 
multigraphs. Connections between the resistance distance and the forest distances 
are established. 
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1 Introduction 



Proximity measures for graph vertices and related algebraic indices have a 
wide range of applications. These fall into information transmission, organic 
chemistry, crystallography, integrated circuit design, urban planning, trans- 
port networks, social networks, politology, aggregation of preferences, epidemi- 
ology, etc. (some references are given in [1]). Let us also mention a few more 
fields of application, namely, cluster analysis (see, e.g., [2]), the theory of par- 
allel computations [3], and optics [4]. In all these areas, there is a need in 
proximity measures different from the classical geodesic distance [5]. 
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Let G be a weighted multigraph with vertex set V(G) = {1, . . . , n} and edge 
set E(G). Multiple edges are allowed, but loops are not; is the weight of 
the pth edge between i and j; the weights of all edges are strictly positive. 2 

In [6] we proposed the relative forest accessibility measure of vertex proximity. 
More exactly, this is a one-parametric family of indices: if a parameter a > 
is chosen, then the matrix Q a = (q^) of vertex comparative proximities is 
given by 

Q a = (I + aL)-\ (1) 



where I is the identity matrix, L = L(G) = (i^) is the Laplacian matrix of 
the weighted multigraph G (also termed as the Kirchhoff and the admittance 
matrix) : 



% = ~J2 £ ij> J^i, h 3 = 1, • • • ,n, (2) 
P =i 

in = £ij, i = 1, ... ,n, (3) 

and is the number of edges with terminal vertices % and j. The matrices 
Q a are doubly stochastic, i.e., 



«5 ^ °' ( 4 ) 

n n 

E ^fe = E 9fcj = = !. • • • > ( 5 ) 

fc=i fe=i 

and symmetric. The entry qf* can be interpreted as the relative share of the 
connections between % and j in the totality of all i's connections with the ver- 
tices of G. The parameter a determines the proportion of taking into account 
long and short routes between vertices. 

In [6] we considered a metric defined by means of the relative forest acces- 
sibilities. The distance between vertices % and j of a (weighted) multigraph 
was defined as q\ + — qjj — q^. In the present paper, we study two para- 
metric families of forest metrics of a graph. The elements of these families are 
proportional, but the properties of the families, taken as a whole, are different. 

Definition 1 For a given parameter a > 0, the value 

d? j = \(Qii + Qj j -Qi j -Qjd, i,j = l,...,n (6) 

2 In this paper, we consider only such multigraphs; sometimes we will simply call 
them graphs. 
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will be called the forest distance between % and j; the value 



Pij = a (<lii + Qjj ~ Q?j ~ Qji), i,j = l,...,n (7) 

will be called the adjusted forest distance between i and j . 

The fact that these functions are metrics follows from Proposition 1 in [6]. 

The multiplier 1/2 in (6) ensures the equality of scales for the relative forest 
accessibilities and the forest distances (see (22) and the subsequent remark). 
If the vertices % and j belong to the same component of G, then Corollary 1 
of [1] implies lim^oo dfj = 0. The asymptotic behavior of the adjusted forest 
distances is studied below in Section 5. In Section 2, we demonstrate that 
the reversed pfj can be used to measure the cumulative weight of connections 
between i and j. 

One of the aims of this paper is to present and interpret formulas for the 
increments of forest distances and relative forest accessibilities under the basic 
multigraph transformations. These expressions are very simple and testify 
to the intrinsic nature of the forest metrics. Also, we will give a stochastic 
interpretation of the forest metric and establish some connections between 
the forest and resistance metrics. 

Until the beginning of Section 5, a will be a fixed parameter, so a simplified 
notation such as d^, p^, q^, and Q will be used instead of the complete 
notation, dfj, pfj, qfj, and Q a , respectively. 

Some properties of the forest metric of a graph have been studied in [6] and [7]. 
It follows from item 5 of Proposition 7 in [6] that if the total weight of all 
edges that have terminal vertices k and t gets an increment of Ae kt > 0, then, 
provided that there are no other changes, the increments of the adjusted forest 
distances in G are 



(Pik -Pit +Pjt -Pjk) 
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Observe that for fixed p ik —p it +Pj t ~Pjk an d Pkti ^ ne increments Ap 4 • are in- 
dependent of a. Besides that, (8) implies that the forest distances cannot 
increase as a result of the addition of a new edge. 

The double stochastic property of Q yields the inequalities 



dij < 1, (9) 
Pij < 2a, i,j = l,...,n, (10) 
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and, by proposition 5 from [6], the equalities in (9) and (10) are attained if 
and only if i and j are two isolated vertices. 



Corollary 9 of [8] provides a more accurate upper bound for gL: 



du < (1 + aa(G)) 



-i 



(11) 



where a(G) is Fiedler's algebraic connectivity of the graph (which is the second 
minimal eigenvalue of L). Then (l+aa(G ! ))~ 1 is, obviously, the second maximal 
eigenvalue of Q. Upper bounds for the algebraic connectivity of the graph can 
be obtained by transforming (11) (see (14) in [7]). 

The diameter of a graph is the greatest distance between its vertices. This way, 
(9) and (11) provide upper bounds for the forest diameter of a graph. Other 
bounds for the forest distances will be given below (equations (16) and (17)). 



2 How forest distances change when connections between vertices 
strengthen 

An expression for the increments of relative forest accessibilities under basic 
graph transformations was given in [6] , item 1 of Proposition 7 (see also related 
results in [7-9]). Now we formulate this property in a more general form. 

Definition 2 We say that a weighted multigraph G' differs from G in a (k, t) 
edge only if for some Ae fci = e ^ 0, G' can be obtained from G by increasing 
the weight of some edge e^ t by e, or by adding a new edge between k and t 
with weight e (e > 0), or by removing an edge between k and t with weight —e 



While saying that G' differs from G in one edge only, by Ae kt = e we will 
denote the weight (or the increment of weight) in the above definition; primed 
notation will relate G', nonprimed notation will relate G. 

Proposition 3 Suppose that weighted multigraph G' differs from G in a (k, t) 
edge only. Then for every i,j e V(G), 



The results of this section rely on the following lemma, which follows from (8) 
after substituting i = k and j = t. 



(e<0). 




Q(gjfc -Qit)(Qjt-Qjk) 
Pkt +£~ L 



(12) 
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Lemma 4 Let G' differ from G in a (k, t) edge only. Then 



11a 

-^r- — = A£ kf (13) 

Pkt Pkt 



Proof of Proposition 3 is carried out by the same argument as the proof 
of item 1 of Proposition 7 in [6]. The only necessary addition is that the 
denominator on the right-hand side of (12) cannot vanish even with a negative 
Ae kt . Indeed, by Lemma 4, 1/ p kt +Ae kt = l/p' kt . As was demonstrated in [6], 
p' kt is nonzero (and exists!) at k ^ t for every multigraph G, consequently, 
VP fc t +^kt 0, therefore, p kt +1/Ae kt ^ 0. 



The connection of an inverse form between the distance and the increment of 
weight in (13) requires some explanation. Note that since Ae kt can be arbi- 
trarily large, whereas the distance between k and t should decrease with any 
increase of Ae kt , a linear connection between the distance and the Ae kt cannot 
be expected. Taking into account that the equality p' kt = p kt is inevitable when 
Ae kt = and that additivity of distance increments and reciprocity principle 
(see below) should also be satisfied, it is difficult to expect for the distance in- 
crements a simpler form than (13). It should be also remarked that the inverse 
connection between the weight of an edge and its length is the most natural 
way of extending the geodesic distance to weighted graphs (cf. [10]). 

The distance (7) for each pair of vertices generally depends on the graph as a 
whole, but if p kt is fixed, then, by Lemma 4, p' kt is determined by Ae kt only. 
In particular, for any unweighted multigraph, the addition of a new (k, t) edge 
yields 

1 1 

v - IT = x > (14) 

Pkt Pkt 

which implies 

Ap kt = -p kt Pkt- (15) 



The form of (13) and (14) suggests the consideration of inverted distances. 

Definition 5 Let i and j ^ i be distinct vertices of G. The value 9^ = 
(Pij)' 1 ~~ (2a)" 1 will be called the cumulative weight of connections between i 
and j in G. 

By Lemma 4, 6' kt = 9 kt +Ae kt , and 9' kt = 9 kt +1 for unweighted graphs. Suc- 
cessively applying Lemma 4 to all the edges with terminal vertices i and j and 
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using (10) along with its equality condition, we obtain the following statement, 
which justifies the term "cumulative weight of connections." 

Proposition 6 Suppose that G is an arbitrary weighted multigraph, i and j 
are distinct vertices of G, and is the total weight of all edges with terminal 
vertices i and j . Then 

1. 9^ = 9®j + £j-, where 9^ is the cumulative weight of connections between 
i and j in the multigraph resulting from G by the removal of all edges with 
terminal vertices i and j. 
2- 9 i3 > £y . 

3. 9^ = if and only if i and j are not connected with the other vertices {but 
can be connected to each other). 

Item 2 of Proposition 6 and Definition 5 provide the following upper bounds 
for the forest distances: 



d ij <(l + 2ae ij )-\ (16) 



Pij <(e ij +(2a)- 1 )-\ i,j = l,...,n, (17) 



where is the total weight of the edges with terminal vertices i and j in G. 
The equality condition is the same as in item 3 of Proposition 6. 



3 Comparing the forest distances in two graphs that differ in one 
edge 



In this section, we study how a modification on a (k, t) edge alters the whole 
profiles of forest distances and relative forest accessibilities. The alteration of 
the values 



T i(kt) — dik—d it , (18) 
T»(fct) = 9ifc-9it. i = l,---,n (19) 

will be of major interest here. 
In particular, we have 

T t(kt) = d kt = —T~k(kt) ■ 

The simultaneous consideration of forest distances and relative forest accessi- 
bilities rises a question: what is the essential difference between them, except 
for the fact that a shorter distance usually corresponds to a greater accessibil- 
ity? The main difference is in the very relative nature of the forest accessibility 
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measure: as follows from Proposition 7 in [6], no distance can increase after the 
addition of a new edge, whereas the increments of a relative forest accessibility 
can have either sign. More specifically, if k is "more accessible" from % than t 
(q ik > q it ), whereas t is "more accessible" from j than k (q jt > q jk ), then the 
addition of a new edge between k and t increases the accessibility of j from i 
(q'ij > because the new edge "widens a road" between % and j. Inversely, 
if from both % and j, k is "more accessible" than t (q ik > q it , q^ k > g-J, then 
the addition of a (k, t) edge decreases the relative forest accessibility of j from 
i {q[j < q^). Here, the new edge connects i and j with t and with the vertices 
situated "beyond" t more intensively than this edge connects % with j (recall 
that the relative forest accessibilities of all vertices from a given vertex sums 
to 1). The diagonal entries of Q measure the "solitariness" of the vertices 
(which is substantiated by Theorem 3 in [7]); no one of them can increase 
when a new edge is added. 

There exists a one-to-one correspondence and a specific duality [11] between 
metrics and S-proximities (one of which is relative forest accessibility). 

Definitions (6), (18), and (19) imply simple connections between d kt , 7i^ kt ^, 
and r m : 



2d kt 


= n k(kt) ~ n t(kt)i 


(20) 


2r i(kt) 


= ( n k(kt) ~ n i(kt)) + ( n t(kt) _7r i(H))' 


(21) 


T i(kt) ~~ T j(kt) 


= n j(kt) ~ n i(kt)i h h k, t — 1, . . . , U. 


(22) 



Remark 7 Equation (22) can be interpreted as the equality of scales of rela- 
tive forest accessibilities q^ and forest distances d^ (see definitions (18) and 
(19)). This equality of scales was attained by putting the coefficient 1/2 in (6). 

The following statement provides simple expressions for the increments of all 
forest distances and accessibilities when one edge is altered. 

Proposition 8 Suppose that G' differs from G in a (k, t) edge only, AcL- = 
d'ij ~~ ^iji an d ^Qij = Qij ~ Qij- Then the following equalities are true for all 
i,j = l,...,n : 



^Qij = -a^i(kt) ^j(kt) -t 1 = - ae A{kt)K){ktyTri ( 23 ) 

a kt a kt 

2Ad. . = -ae(r m) -r m) ) 2 ^ = -ae{r' l{kt) - r' j{kt) ) 2 ^ 

a kt a kt 
d! d 
= -<xe(ir i{kt) = -»« { kt) ~ ^W 2 /' ( 24 ) 

a kt a kt 

where e = Ae kt is the weight of a new edge or the increment of the weight of 
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an existing edge between k and t. 



The first equalities in (23) and (24) are proved by the application of Lemma 4 
to (12) and (8), respectively. The third equality in (24) follows from (22). 
The other connections can be derived by applying the following reciprocity 
principle: 



Reciprocity principle. Suppose that a statement A is true for every pair of 
multigraphs (G, G') such that G' differs from G in a {k, t) edge only. Suppose 
that statement A' results from A by replacing all values relating to G with 
the corresponding values relating to G' and vice versa and interchanging Ae kt 
with Ae tk . Then A' is also true for every pair of multigraphs (G, G') such that 
G' differs from G in a (k, t) edge only. 



This principle resembles the oriented duality principle, but its nature is differ- 
ent. The reciprocity principle follows immediately from the following obvious 
fact: if G' differs from G in a (k,t) edge only, then G also differs from G' in 
a (k,t) edge only, but with the reversed increment of weight. The reciprocity 
principle, in spite of its obviousness, can be used to obtain some nonobvious 
corollaries (for example, the equalities in Proposition 8). 

Let us specify now how the values T^ kt ^, n^ kt ) change under the basic trans- 
formations of a multigraph. 

Proposition 9 Let weighted multigraph G' differ from G in a (k, t) edge 
only. Then for all i = 1, . . . , n, the following equations hold: 



T i(kt) _ T i(kt) 

d'kt d kt 
d'kt d kt 



or, equivalently, ik u = ik u ; (25) 
d-kt d kt 

or, eqmvalently, q ' ik ~ — = g<fc ~ git . (26) 



Equation (26) is a corollary from Proposition 3 and Lemma 4, but it can be 
alternatively obtained from the equations in (23) and (24) that are proved by 
means of reciprocity principle; (25) follows, say, from (26) and (21). 

According to (25), the difference of distances d ik —d it reduces as k and t come 
closer, 3 so as the proportion in (25) is preserved, which makes sense; according 
to (26), the same property holds true for the difference of accessibilities. 

3 This "coming closer" is caused by strengthening connections between k and t, 
whereas the other connections remain the same. 
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Observe also that, by (25)-(26) and Lemma 4, the ratio 



T 



i(kt) 



i(kt) 



TV 



i(kt) 



i(kt) 



4 

d. 



kt 



kt 



pk = 1 
pkt 1 + £ Pm 



1 - zp'kt 



(27) 



is the same for alH = 1, . . . , n. 



4 An interpretation of the forest metric 

In this section, we give an interpretation of the forest distance between % and 
j in terms of spanning rooted forests that connect i and j. 

Recall that a rooted forest is an acyclic graph with one vertex marked as a root 
in each its component. The components of a rooted forest are rooted trees. 

Every spanning rooted forest of G where % is a root of a tree and j does not 
belong to this tree will be called an unsuccessful connection from i to j. A 
successful connection from i to j is a spanning rooted forest of G where j 
belongs to a tree rooted at i. 

By G a we denote the weighted multigraph obtained from G by multiplying 
the weights of its edges by a. 

Consider the following model of choosing a random connection. 

A stochastic model of connecting vertices % and j. 

1. Choose one of the two vertices, % or j, with probability 1/2 each. 

2. Choose a spanning rooted forest of G a at random: the probability of choos- 
ing forest F is 



where jF(G a ) is the set of all spanning rooted forests of G a and e(F') is the 
weight of the forest F', which is defined as the product of the weights of its 
edges. 

3. If the forest chosen at stage 2 is an unsuccessful connection from the vertex 
chosen at stage 1 to the remaining vertex of the pair then we say that 

an unsuccessful connection between i and j is chosen. 



p(F) 



E e(F'Y 



(28) 
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The following statement is an interpretation of the forest distance between i 
and j. 

Proposition 10 The forest distance is equal to the probability of choosing 
an unsuccessful connection between i and j in the above model of connecting i 
and j . 

Proposition 10 clarifies the concept of forest metric: two vertices are close to 
each other if the probability of choosing an unsuccessful connection between 
them is small; the distance is precisely equal to this probability. 

PROOF. Observe that, by the matrix-forest theorem [6], — g^- is the 
weighted fraction of unsuccessful connections from i to j. Therefore, by (28), 
Qu 1S the probability of choosing an unsuccessful connection from i to 
j on stage 2 of the model. Consequently, by the total probability formula, 
\{o.u + liQjj i s t ne probability of choosing an unsuccessful connec- 
tion between i and j, and, by (6), this value is equal to d^. 

Making use of Theorem 3 in [1], an interpretation of this kind can be for- 
mulated for the resistance distance too. Connections between the resistance 
metric and the forest metrics will be considered in the following section. 



5 Relations between the resistance metric and the forest metrics 

According to item 3 of Proposition 6, the index Q { - referred to as a cumulative 
weight of connections between i and j has a rather strange property: if i and 
j belong to different components of G, but at least one of these vertices is not 
isolated, then 9^ > 0. Which connections can be associated with such a pair 
of vertices? The following analysis of relations between the resistance metric 
and the forest metrics clarifies the point. 

The resistance metric can be defined as the function that associates a nonneg- 
ative number on the extended real line 1RU {+00} with every pair of vertices 
i, j G V{G): if i and j belong to different components of G, then p { - = +00; 
otherwise 

p,i (;,-(jj-(u-(ji- (29) 

where £+• ,i,j = l,...,n, are the entries of the Moore-Penrose generalized 
inverse L + of the Laplacian matrix L. An electrical interpretation and some 
properties of the resistance metric can be found in [12-15]. 



10 



Consider the asymptotic behavior of the forest metrics d a and p a . As a — > 0, 
the metric d a tends to the discrete metric: 

«& = {!: i;t (3o> 



whereas p a tends to the zero function. 

Let V i be the vertex set of the component of G that contains vertex i. As has 
been mentioned in Section 1, if j e V \, then df 3 ; — > as a — > oo. The behavior 
of o? Q and p a as a — > oo is completely specified by the following statement, 
which follows from Corollary 1 and Proposition 9 of [1]. 

Proposition 11 The limits d°° = lim^oo dfj and p°° = lim^^ pf^ always 
exist (in the latter case, on the extended real line), and 



d°° 



?Ti = P.. 



0, | ( jev t , 



iji 



i,j = l,...,n. 



(31) 
(32) 



Note that if j £ V { , then the limit of the cumulative weight of connections 8^ 
as a — > oo is 0. 

Another relationship between the forest metrics and the resistance metric is 
also interesting. 

Definition 12 A weighted multigraph G' will be called the a-extension of the 
weighted multigraph G if 

(i) V(G') = V(G) U {0} ; and 

(ii) the restriction of E(G') to V(G) coincides with E(G), and the ratio of the 
corresponding edge weights in G' and in G is a, and 

(Hi) E(G') contains one edge (0,i) of weight 1 for each vertex i e V{G). 

Proposition 13 For every weighted multigraph G and every i,j G V(G), 
= 2d?j(G) is true, where G' is the a- extension of G, p^(G') is the 
resistance distance on G' , and dfj(G) is the forest distance on G with param- 
eter a. 



PROOF. 1. The matrix Q Q a obtained from Q a = (I + aL)^ 1 by the addition 
of the zero row and zero column corresponding to the vertex is a generalized 
inverse of L' , the Laplacian matrix of G' . Indeed, the required equality 

L'Q° a L' = L>, (33) 
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is verified straightforwardly by the multiplication of block matrices. 

2. As has been demonstrated in [16] (see also [17]), every generalized inverse 
H of a doubly centered matrix Y with singularity 1 can be represented as 



where Y + is the Moore-Penrose generalized inverse of Y, a and b are some 
vectors, and e = (1, 1, ... , 1) T . 

3. Observe that the transformation 



performed with any matrix H = (h^) that satisfies (34) and with Y + gives the 
same result. Applying this observation to the matrix L' which has singularity 1 
and making use of (33), we obtain the required statement. Mention that the 
identity of the results of transformation (35) for all generalized inverses of the 
Laplacian matrix of a connected graph has been noticed in [14]. 

We are now in position to discuss the above-mentioned property of the cu- 
mulative weight index 0^. this weight can be positive when % and j belong to 
different components of G; on the other hand, its limit value corresponding to 
the resistance distance is zero. 

The objects modeled by the vertices of the same graph usually have a specific 
similarity or a common derivation, which causes their a priori relationship. The 
forest metrics enable one to take such a relationship into account in the model 
and to control the strength of the a priory connections by the parameter a: an 
increase of a corresponds to a weakening of a priory connections. From this 
point of view, considering the value (p^) -1 , instead of 6^ = (Pij)^ 1 — (2a) _1 , 
as the cumulative weight of connections between % and j makes sense too. In 
accordance with Proposition 13, it is not additional connections between the 
objects that represent their a priori relations, but connections with the intro- 
duced "hidden source" vertex. This way, the input structure of connections 
between the objects remains intact. 

The introduction of a "hidden source" frequently enables one to obtain simpler 
proofs for graph-theoretic statements. For instance, the set of spanning rooted 
forests of any multigraph G can be put into a one-to-one correspondence with 
the set of spanning trees of the 1-extension of G; the corresponding subgraphs 
have the same weight. This method makes it possible to reduce some theorems 
that involve rooted forests to theorems about trees. 



H = Y + + ae T + ea T , 



(34) 




K +hjj -/',, -hji, 



(35) 
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In some models, varying the parameter a of a-extension enables one to specify 
the relative importance of short and long connections in a graph. For instance, 
this is the case for the linear regression analysis of paired comparisons, where 
the addition of a "hidden source" corresponds to the turn from the least 
squares estimates to the ridge estimates of the object effects. The matrix of 
the system of linear equations corresponding to the ridge estimates is non- 
singular as distinct from the matrix of normal equations. In other words, a 
regularization procedure applied to the ill-posed least squares problem leads 
to the ridge estimates for the object effects. It turns out that ridge estimates 
are more adequate in some cases than the least squares estimates of the object 
effects in paired comparisons [18]. These ridge estimates can be represented 
via relative forest accessibilities [19]. 

The simple facts and relationships presented in this paper contribute to the 
analysis of the resistance and forest metrics and provide some framework for 
the comparison of various graph metrics. 
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